Selective focusing of electrons and holes in a graphene-based superconducting lens 
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We show that a graphene pnp junction with a central superconducting electrode acts as a Veselago 
lens for incoming electrons by focusing them and their phase-conjugated counterpart (holes) into 
different points of the optical axis. This selective focusing suggested by a simple trajectory analysis is 
confirmed by fully microscopic calculations. Although the focusing pattern is degraded by deviations 
from the ideal conditions we show that it remains visible for a wide range of parameters. We discuss 
how this property can be useful for the detection of entangled electron pairs. 
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I. INTRODUCTION. 

The possibility of fine-tuning the density of carriers 
in graphene together with their resemblance to massless 
particles like photons have made graphene a promising 
candidate for testing photonic analogies in electron trans- 
port. Owing to the relativistic chiral nature of carriers 
there is a suppression of backscattering at a graphene- 
based pn junction known as Klein tunneling _F, which 
has been confirmed in recent experiments [2J [3J. Fur- 
thermore, an electron beam flowing through a single pn 
junction experiences negative refraction [5] and thus this 
system has been proposed as the electronic equivalent 
of a Veselago lens [5]. Thanks to this analogy, striking 
properties of the meta-materials such as perfect-lensing 
[5] could be explored in graphene. In particular, super- 
collimation of electron beams [7J has been proposed for 
graphene under periodic potentials and the focusing of 
electron beams has been studied in circular graphene 
injunctions jS], graphene nanoribbons [9] and the sur- 
face of topological insulators [TO] . Graphene-based Vese- 
lago lenses have also been proposed as filtering systems 
for spin-polarized electron beams [UJ. Additionally, ad- 
vances in the construction of ballistic pnp junctions in 
graphene have been reported [T2Ml5| . 

New interesting possibilities can emerge if we consider 
the case of hybrid graphene-superconductor nanostruc- 
tures [TBHTS] . Good contact can be achieved between 
lithographically defined superconducting electrodes and 
graphene layers [T5H22] . In such devices a super- 
conducting gap is induced by proximity effect on the 
graphene region underneath the metallic electrodes; in 
these conditions Andreev processes featuring conversion 
of electrons into holes take place [23] • In a graphene- 
based normal-superconductor-normal (GSG) junction, 
local and crossed (CAR) Andreev reflections can occur 
if the width of the central superconducting electrode 
is comparable to the superconducting coherence length 
[23] . The time-reversal of these processes corresponds to 
the splitting of a Cooper pair from the superconductor 
into an entangled electron pair in the normal electrodes 
|25j . Although progress has been achieved in the exper- 
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FIG. 1: (Color online) Graphene superconducting pnp junc- 
tion. The band structure of each region is shown on top for 
the case when the normal electrodes L and R are adjusted 
to have the opposite doping level than the superconducting 
central region S. Incoming electrons from the left electrode 
(solid red lines) transform into evanescent electron and hole- 
like excitations inside the superconductor (dashed red and 
blue lines, respectively). In the right region, electrons and 
holes (solid blue lines) focus into distinct regions. 



imental realization of Cooper pair splitters using carbon 
nanotube and semiconducting nanowire quantum dots 
[2T51 [27] . it is expected that graphene can provide even 
better conditions for the entanglement detection. 

In this letter we propose to create a Veselago lens in 
a GSG junction to focus electrons and holes in differ- 
ent spatial regions. The idea is schematically depicted 
in Fig. [I] A graphene sheet is deposited on top of two 
independent gate electrodes and a central superconduct- 
ing electrode (denoted L, R and S respectively). The 
superconducting electrode shifts the electronic bands of 
the underlying graphene region by transference of elec- 
trons to produce an n-doping effect. Subsequently, the 
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FIG. 2: (Color online) (a) Map of the electron transmission 
T ee (x,y) near the focusing point x = 3Ws/2. The width of 
the superconducting region is Ws — £ ~ 1/ttm and the doping 
levels are |i?^' S,ii | = 0.5eV. (b) For the same parameters, 
map of the hole transmission T e h(x, y) near the focusing point 
x — 5Ws/2. (c) Sketch of the trajectories of electrons and 
holes through the superconductor into the region R. See the 
text for details, (d) Logarithmic scale plot of the transmission 
of electrons (red) and holes (blue) into region R for y — 0, 
showing the profiles of (a) and (b) where the peaks due to 
internal reflections can be distinguished. 



doping level of each normal region is adjusted by the gate 
electrodes so the system behaves as a pnp junction but 
with the peculiarity that superconductivity is induced in 
the central region. Injection of electrons can be real- 
ized in the L region by means of a local probe and when 
the central electrode is in the normal state the system 
acts as a Veselago lens focusing electrons at the region 
R. When superconductivity is "switched on" , evanescent 
electron and hole-like states are created in the graphene 
region located under the central electrode. In the right 
electrode, electrons and holes become propagating waves 
again and are focused into regions separated by hundreds 
of nanometers. This spatial separation would allow the 
detection of the transmitted holes by means of a second 
local probe (detector in Fig. [I]). 



II. INDEPENDENT FOCUSING OF 
ELECTRONS AND HOLES. 

For modeling the Veselago lens depicted in Fig. [I] we 
consider an infinite plane of graphene with a supercon- 
ducting electrode covering the region < x < Ws, while 
the regions x < (L) and x > Ws (R) remain in the 
normal state. Although graphene is not intrinsically su- 
perconducting, the superconducting electrode can induce 
a pairing amplitude A by proximity effect [28-30] . We 



choose the width of the central region to be comparable to 
the superconducting coherence length Ws — £ = hvp/A. 
For a typical superconductor like Pb or Al, A ~ lmeV 
and thus £ ~ 0.5 — 1/zm. An Al or Pd/Al electrode in 
the superconducting state induces a n-doping in the un- 
derlying graphene region estimated as Eg ~ — 0.5eV [31]. 
Since the gate potentials of the normal regions can be ad- 
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justed independently we choose them to be E F = —E%. 
This set of parameters define a pnp junction with a cen- 
tral superconducting electrode. The transport properties 
between the two normal electrodes are computed in terms 
of one-particle Green functions following the method ex- 
plained in Ref. |32j and briefly presented in Appendix 
[A] This allows to define an electron transmission proba- 
bility T ee (x,y) from region L to region R (i.e. electron 
cotunncling EC), and a CAR probability T e h(x, y) where 
[x, y) denotes the detector coordinates when the source 
coordinates are (— d, 0). In Fig. [2] we show the result 
of the microscopic calculation for the transmission prob- 
abilities. T ee exhibits a well-defined peak at 3<i while 
we obtain a maximum of T e h at 5c?, as it is shown in 
Fig. [2ja) and Fig. (2jb) respectively. Both results are 
normalized to To = T ee (—d, 0). The intensity of the elec- 
tron focusing is much higher than that of the holes, i.e. 
T ee (3o!,0) - 100T e/l (5d,0), thus indicating that the EC 
signal is much greater than the CAR signal. However, the 
spatial separation between peaks is exactly the width of 
the superconducting region, which is of the order 1/im. 
This would allow to detect each signal independently as 
we discuss in more detail below. 

These results can be qualitatively explained with a sim- 
ple analysis of the group velocities of the particles at each 
region (see Fig. [2][c)). For a perfectly symmetric pnp 
junction, the normal regions L and R are p-doped with 
Ep > 0, while the superconducting region S is n-doped 
with —Ep. In the heavily doped regime \Ep\ 3> A,E 
an analysis of the propagation of waves based on clas- 
sical trajectories is sensible since Xf <C Ws- Even for 
the evanescent waves within the superconducting region 
and with \E\ < A this type of analysis is valid consid- 



ering waves with kp ~ k 



L,R 



Ep/hvp. Subsequently, 



the angle of incidence and transmission of particles is de- 
fined as (/>(k y , Ep) = ± a,TCsm[hvpk y /Ep], where the sign 
depends of the doping level of each region. As a result, 
the group velocities of electrons and holes can be written 
as V ej /j = ±£VFke,h/\^e,h\ = ±£Vf (cos 0, sin 0) where 
e = 1(— 1) for quasiparticles in the conduction (valence) 
band. From the conservation of the component of the 
wave vector parallel to the interface we reach the elec- 
tronic equivalent of Snell's law at each interface |U [TO] 
which allows to define a relative refraction index at each 
interface Uls.sr — sm(f> L ' S / sin(p S ' R = —e. Taking into 
account the sign due to the band index (conduction or 
valence) and the one due to particle index (electron or 
hole), when the particle type is the same at both sides 
of the interface (e = 1), the change of band causes a 
negative refraction with n = — 1. On the other hand, 
when the particle type is not conserved (s = —1), there 
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FIG. 3: (Color online) Microscopic calculation for T ee (x,0) 
and T e h(x,0) when the injection of electrons is done at d = 
— 3Ws/2. The inset shows a sketch of the classical trajecto- 
ries. Notice that peak of T e h at x — 3.5Ws is ~ 5 times larger 
than the EC background. 



is no negative refraction and n = 1. This explains the 
classical trajectories sketched in Fig. [I] and Fig. [2jc) : 
incoming electrons from the region L (solid red lines) 
transform into electron and hole-like excitations inside 
the superconductor (dashed red and blue lines, respec- 
tively). While the former experiences a negative refrac- 
tion because the particle type is conserved, the latter 
follows the same path as the incoming electron. At the 
second interface, the processes that preserve the particle 
type experience a negative refraction and are focused to 
form an image in the optical axis x (for an analysis of 
the transmission amplitudes see Appendix |Tj| . 

Furthermore, at each interface specular reflection oc- 
curs when the particle type is conserved while retro- 
reflection happens otherwise (Fig. [2jc)) . Electron and 
hole-like excitations can endure two consecutive specu- 
lar reflections inside the superconductor to create a new 
electron or hole beam in the normal region R. This leads 
to a sequence of alternated electron and hole focusing 
points at the optical axis. In Fig. [2jd) we show the mi- 
croscopic calculation in which we obtain peaks of T ee at 
x = (2m — \)W$ and peaks of T eh at x = (2m + \)W~s, 
with m = 1,2,.... The intensity of these peaks decays 
exponentially with the distance to the superconductor, 
consistently with the behavior of the proximity effect in 
a graphene-superconductor interface [2"BH3TI] . 

While the previous analysis explains the separate fo- 
cusing of electrons and holes, it would be desirable that 
the background EC conductance at the CAR peak could 
be further reduced. This can be achieved by increas- 
ing the injection distance d. Indeed, the injection point 
(— d, 0) determines the origin of the sequence of focus- 
ing points in which the separation between maxima of 
T ee and that of T e h is Ws- For d < Ws the sequence 




FIG. 4: (Color online) (a): Map of T e h(x,y) for an asym- 
metric pnp junction. We set Ep = 0.4eV and keep the rest 
of the parameters the same as in Fig. [3] where d > Ws- The 
refracted waves in the right region form a characteristic inter- 
ference pattern showing two caustic curves near the cusp for 
each of the transmitted holes, (b): Map of T e h(x, y) for a pnp 
junction calculated using the TB model. The length of the 
central superconducting region is Ws = 85.2nm. The gate 
potentials used are Ep' R = —Ep = 0.6eV, with A = 2.7meV. 
The left panel includes no disorder. The central panel has 
Vo = 54meV. The right panel has Vo = 81meV. The smearing 
of the potential for the three panels has a range of ~ 35nm. 



of points starts always with a maximum of T ee . On the 
other hand, for d > Ws, electrons are only transmitted 
into the region R after two or more internal specular re- 
flections. However, a combination of one specular reflec- 
tion and one retro-reflection, which changes the particle 
type, allows to have a hole focusing point at x = 1.5Ws, 
i.e. to the left of the electron focusing point at x = 2.5Ws 
(see Fig. [3] for a sketch of the trajectories and the mi- 
croscopic calculation). The intensity of this peak is the 
same as the CAR peak appearing at x = 3.5Ws, to the 
right of the electron focusing point, which comes from a 
transmitted hole without internal reflections. Both peaks 
of T e h in Fig. [3] are greater than the peak of T ee and in 
particular the peak at x = 3.5W$ is ~ 5 times larger 
than the background EC contribution. 
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III. DEVIATIONS FROM THE IDEAL CASE. 

We have considered thus far that each interface is a 
perfectly symmetric injunction, i.e. E P F = —E F . If the 
doping level of one of the regions is not perfectly aligned 
with the next we have that for that interface. In 

Fig. |Z|a) we show a map of T e h(x, y) with E F = —Ep = 
0.5eV and E§ — 0.4eV where the refraction index of 
the SR interface becomes n = —Ep/Ep = —0.8. For 
this case the classical trajectories are deformed at the 
region R. The focal points are displaced from the ones 
shown in Fig. [3] and the envelope of the refracted rays 
becomes a caustic curve |33j . The focusing pattern is 
reproduced but the intensity of the cusps is not the same, 
as it was the case for the symmetric junction. Caustic 
curves for electrons in asymmetric injunctions have been 
predicted to appear [U [51 [TU] but we show here that these 
curves appear in spite of being originated from hole-like 
excitations inside the superconductor. We thus conclude 
that a doping imbalance between regions has the same 
effect on the CAR signal that on the EC one. 

The results presented thus far correspond to an infinite 
layer of pristine graphene under a perfectly sharp poten- 
tial profile. A more realistic model should include size 
effects such as a graphene layer with a finite length, a po- 
tential profile varying smoothly along the sample and the 
inclusion of electron-density inhomogcneitics (i.e charge 
puddles [34l [35])- Using a tight-binding (TB) model we 
explore the stability of the focusing pattern under these 
premises. We define a defect-free graphene strip of to- 
tal length W ~ 900nm, with a central superconduct- 
ing region of length Ws ~ 90nm and coupled to normal 
metallic electrodes at the edges (the details of the TB 
model are presented in the Appendix [C| . A smearing 
of the potential profile within a range of 30 — 40nm is 
introduced at each normal-superconductor interface. In 
addition, we introduce random inhomogeneities of the 
potential profile of strength Vq over an area of typical 
length d ~ 20 - 30nm. 

We show in the left panel of Fig. [5|b) the TB re- 
sults for T e h(x,y) in the absence of disorder (Vq = 0). 
The focusing spot at x = 2.5Ws is clearly distinguish- 
able although some diffraction effects are present due 
to the finite length of the system and the smearing of 
the potential. When we introduce disorder of strength 
Vq = 54meV (central panel) the intensity of the focusing 
spot is reduced but its size remains almost unchanged. In 
the right panel, when the disorder strength is increased 
to Vo ~ 81meV, diffraction effects overcome the focusing 
pattern. The smearing of the potential has a range of 
~ 35nm for the three panels. As it was demonstrated 
in Ref. [9] for a smooth pn junction, the smearing of 
the potential introduces a small diffraction effect in the 
transmission, reducing the intensity of the focusing point, 
but leaving the extension of the spot almost unchanged. 
We also find that the focusing pattern is robust against 
disorder caused by charge puddles of width ~ 20 — 30nm 
and strength Vq < 80meV. These parameters are well 



above the measured inhomogeneities in graphene, which 
are bounded to 30meV over 20 — 30nm j34] [35] . 



IV. CONCLUSIONS. 

In conclusion we have shown that selective focusing 
of electrons and holes can be produced in a GSG junc- 
tion. In addition, the geometry can be tuned in or- 
der that the CAR peak dominates over the EC back- 
ground. We have also shown that the focusing is robust 
against deviations from the ideal conditions. Under these 
premises, a possible experimental realization of this pro- 
posal is sketched in Fig. [I] where electrons are injected 
at the source electrode and collected at the detector elec- 
trode. If the source is a fixed electrode, a mobile detector 
would be able to distinguish between the EC signal and 
the CAR signal by moving from one focusing point to 
the other. On the other hand, if the detector is fixed, 
by moving the source electrode the focusing pattern can 
be adjusted to reach the fixed electrode. Although the 
nonlocal electron-hole transmission is reduced by a fac- 
tor 10 -4 with respect to the transmission at the injec- 
tion point, the total nonlocal conductance can reach a 
measurable value when adding the contribution of many 
channels. These properties open an interesting route for 
the detection of entangled electron pairs over distances 
of ~ 1/im. 
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ported by COLCIENCIAS, project 110152128235 (SG 
and WJH) and MICINN-Spain via grant FIS2008-04209 
and EU project SE2ND (PB and ALY). 



Appendix A: Modeling the system. 

We consider an impurity-free graphene sheet in the x — 
y plane. A superconducting electrode is deposited on top 
of the region < x < Ws- The low-energy excitations 
of the system are described by the Dirac-Bogoliubov-de 
Gennes (DBdG) equations 

(H-V(x) A(x) \fu\_ f u\ ,. n 
^ A(x) V{x)-H J\v ) \v )' (AL) 

where H = Hvp {a x k + <J y q) is the one particle Dirac 
Hamiltonian with Fermi velocity vf, A(x) is the pair- 
ing amplitude, V(x) is the potential profile and E > 
is the excitation energy. We impose rigid boundary 
conditions at the normal-superconducting interfaces to 
the pairing and the electrostatic potentials such that 
A(\x\ > Ws) = 0, A(\x\ < W s ) = A and V(x < 
0) = E%, V(0 < x < W s ) = E F , V(x > Ws) = 
Ep. Whenever the pairing potential is assumed con- 
stant and non-zero, the low energy spectrum is given 

by E = J A 2 + (Ep - hv F y/k 2 + q 2 Y. We define the 
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transversal momentum as hvpk± 



with Q — v 1 E 2 — A 2 and Hq the conserved momentum 
parallel to the interfaces. The pairing potential couples 
electrons and holes from different valleys. Eq. |A1| is 
therefore written in Nambu and pseudospin space, omit- 
ting the valley and spin degeneracies. 

The transport properties can be expressed in 
terms of one-particle Green functions which satisfy 
[(E ± i0+) - H(x)] G r ' a {x, x') = S(x - x'), where H de- 
notes the full Hamiltonian of the left hand side of Eq. 
Al| We calculate the Green functions by solving sep- 
arately each region and combining the results following 
the method explained in Ref. [32]. To fully resolve spa- 
tially the Green functions we use the Fourier transform 
g(x,x';y-y') = J dqe iq ( y - y ">G(x, x'). Therefore, by set- 
ting the electron injection in the left region (V < 0) we 
analyze the electron transmission into the right region 
defining T ee oc Tr \Q ee [x, x'\ y — y')\ 2 and the CAR prob- 
ability as T e h Tr \Q e h(x, x'; y — y')\ 2 , where the trace is 
done in the pseudospin space. 



Appendix B: Scattering amplitudes at the 
graphene-superconductor interface. 

We consider a normal-superconductor interface along 
the y-direction in a graphcne sheet, with the normal re- 
gion extended at x < 0. An incoming electron into the 
interface from the normal region can be reflected as an 
electron or a hole with probability amplitudes r ee and r e h 
respectively or it can be transmitted into the supercon- 
ducting region as an electron-like or a hole-like quasipar- 
ticle, with probability amplitudes t ee and t e h respectively. 
The scattering states in both regions are 



+r eh e ik * x 
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where the dependence on y has been omitted because 
the vertical momentum q is conserved. Following the 
notation explained in Ref. [32) . we have defined the 
bispinors in sublattice space ipf e lh = (1, e e tp^ e 2h = 
(1, — e =F4Qc fe ) and the BCS coherence factors u 2 (v 2 ) = 
(E±Q)/2A, which are normalized so that |u| 2 + |v| 2 = 1. 
We are interested in the heavily doped regime with 
\Ep\ = \Ep\ = Ep A,E, which satisfies the mean- 
field approach for superconductivity. In this regime 
« kf h and k^- s = k^' S . The angles at each region 
are thus defined as — = <j) and af — a? — s<p with 
cf> = arcsin Hvf/Ef and s = sign(Ep)sign(Ep), which is 
positive when there is no change in the doping level at 
the interface and negative otherwise (i.e. for a np or a pn 
junction). Matching the scattering states at the interface 



it is straightforward to obtain the reflection and trans- 
mission amplitudes (a detailed description is given in Ref. 
[IB]), in the energy regime of this work they reduce to 
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(B5) 
(B6) 



When the NS interface is equivalent to a perfectly trans- 
parent symmetric nn junction the incoming electron is 
transmitted into the superconductor only as an electron- 
like quasiparticle since s = 1 and thus = 0. 

On the other hand, we can analogously examine the 
reflection and transmission amplitudes of an incoming 
electron-like excitation from the superconductor into the 
rightmost normal region. The scattering states are thus 
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The resulting amplitudes, in the energy regime used are 
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The incoming electron-like excitation can be reflected in- 
side the superconductor both preserving and changing 
the particle type. The former is a specular reflection with 
a change of sign in both components of the velocity while 
the latter is a retro-reflection in which the new excitation 
follows back the path of the incident one. Equivalent re- 
sults are obtained for an incoming hole-like excitation. 
As in the previous case, for a injunction (s = 1) there is 
only transmission into the normal region preserving the 
particle type since t' eh = 0. 

To summarize, for a pnp junction the incident electron 
from the normal region L splits into electron and hole-like 
excitations inside the superconductor and this splitting 
determines the focusing points of each type on the nor- 
mal region R. On the other hand, for a nnp junction 
only electron-like excitations are created inside the su- 
perconductor and thus there is only focusing of electrons 
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in region R. This condition can be relaxed if the trans- 
parency of the interface is not perfect or if the junction 
is not symmetric (i.e. \Ej?\ ^ 



\E S F \ 



Appendix C: Tight-binding model. 

In order to describe a more realistic system, we analyze 
the electronic states of a defect-free graphene layer using 
the tight-binding approximation, 



H = -t, 



E 

<ij> 



(CI) 



where t g = 2Hvf/3(Io ps 2.6eV denotes the hopping ele- 
ment between nearest carbon atoms on the hexagonal lat- 
tice, cio ~ 0.14nm is the smallest carbon-carbon distance 
and Vi is the potential applied to the lattice. The spin de- 
gree of freedom has been omitted due to degeneracy. We 
assume a finite horizontal length W — Wp + W$ + Wr, 
where Wl,r are the lengths of the normal regions and 
Ws is the length of the superconducting one. We have 
well-defined zigzag edges along the vertical direction. We 
impose periodic boundary conditions in the direction par- 
allel to the edges in order to work in the regime in which 
the vertical distance is much larger than the horizontal 
one and thus have a continuum of transversal modes. 



To compute numerically the transport properties we 
connect the zigzag edges to heavily doped graphene leads, 
maintaining the graphene sublattice structure at the 
edges and thus representing the experimental situation in 
which the electrodes are deposited on top of the graphene 
layer [29, 36 40j. The self-energies on the graphene sites 



at the layer edges are thus given by Fp, — Fr — t„. Anal- 
ogously, the central superconducting region has an effec- 
tive coupling with the superconducting electrode given 
by T s - A - lmeV. 

The electrostatic potential profile along the graphene 
strip is defined as 



V(x) = E S F + 



Ep — Ep 



W L 



E§-E S F 



-arctan 



arctan — 

x-W L -Ws 
aa 



(C2) 



L Ft 

where E F ' are the gate potentials of the normal regions, 
E§, is the doping level of the superconducting region and 
the parameter a controls the smearing of the potential at 
the interfaces. In addition, we introduce an extra term 
8V in the potential. This term accounts for disorder in 
the distribution of charge on the graphene strip. SV takes 
random values in the range [— Vb, Vb] over an area of typ- 
ical size d. 
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